Long-range Ordering of Vibrated Polar Disks 



o 

Oh 
< 



o 

s 



o 



> 

o 

o 

CO 



^ 



C. A. Weber/ T. Hanke/ J. Deseigne,^ S. Leonard,^ O. Dauchot,'^'^ E. Frey,i and H. Chate^'^ 

^Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, 

Department of Physics, Ludwig- Maximilians- Universitdt Milnchen, 80333 Munich, Cermany 

'^ Laboratoire de Physique, ENS de Lyon UMR 5672, 69007 Lyon, France 

^Service de Physique de VEtat Condense, CEA-Saclay, URA 2464 CNRS, 91191 Gif-sur-Yvette, France 

^EC2M-Culliver, ESPCLParisTech and CNRS UMR 7083, 75005 Paris, France 
^ Max Planck Institute for the Physics of Complex Systems, Nothnitzer Str. 38, 01187 Dresden, Cermany 

Vibrated polar disks have been used experimentally to investigate collective motion of driven 
particles, where fully-ordered asymptotic regimes could not be reached. Here we present a model 
reproducing quantitatively the single, binary and collective properties of this granular system. Using 
system sizes not accessible in the laboratory, we show in silico that true long-range order is possible 
in the experimental system. Exploring the model's parameter space, we find a phase diagram 
qualitatively different from that of dilute or point-like particle systems. 
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Collective motion in driven or self-propelled particle 
systems is a topic of recent interdisciplinary interest [3]- 
[6]. Within physics, following the works of Vicsek et 
al. [71 [8] and Toner and Tu [9HTT], most progress was 
achieved by studying microscopic models [71 [3 [T2H24] 
and their continuous descriptions J9 HTT| [25H38] . For the 
simplest situation in which the surrounding fluid can 
be neglected ("dry flocking") and the sole interaction is 
some local effective alignment, a picture of basic univer- 
sality classes has emerged, which connects models simi- 
lar to the Vicsek model [7 to continuous theories of the 
Toner-Tu type jOHTT] [25H3T] . Among the landmark re- 
sults are the possibility of true long-range orientational 
order in two dimensions, the generic presence of strong, 
long-range correlations [9 HTH [25] and/or spontaneously 
segregated dense and highly ordered nonlinear structures 
in moving, ordered, fluctuating phases [291 I3T ] . 

These numerical and theoretical results still largely 
lack experimental confirmation. This is mostly due to the 
fact that decisive experimental tests must be performed 
on large numbers of objects under controlled conditions. 
The advent of experiments using purified proteins (mo- 
tors, filaments, etc.) offers a promising playground [39j- 
[43j, but another line of attack, for dry flocking, is to build 
on the experience of the granular physics community, and 
to shake man-made objects p^ [2| I44H48] . Recently, some 
of us have designed and studied the collective motion 
of vibrated polar disks, i.e. millimeter-size objects with 
a built-in oriented axis and a circular top metallic part 
rendering the particles isotropic with respect to collisions 
(Fig. [l^,b; [flu])- Large-scale collective streams and 
anomalous, "giant" number fluctuations were reported 
in collections of approximately a thousand disks mov- 
ing on a carefully vibrated plate. Unfortunately, in this 
experiment — as in others involving man-made objects 
[44-48 — the number of particles used was still too small 
to reach asymptotic results. Moreover, the most ordered 
regimes that could be explored were probably close to 



the onset of collective motion, making it impossible to 
disentangle the properties of the ordered moving phase 
from those of the order-disorder transition. 

In this work, we bypass the inherent difficulties of the 
experimental setup for vibrated polar disks by studying 
the system in silico: We construct a model for the mo- 
tion and collisions of the polar disks of pj |2J , which ac- 
counts quantitatively for most of the known experimental 
properties at the single and pair interaction level. Our 
model also agrees well with observations at the collective 
level. This allows us to study system sizes and bound- 
ary conditions unreachable in the laboratory. We show 
that even in the most ordered regimes observed exper- 
imentally, no long-range collective motion exists. How- 
ever, changing parameters only slightly we find ordered 
regimes which could be observed in the laboratory. Ex- 
ploring the model's parameter space systematically, we 
discover a phase diagram qualitatively different from that 
of dilute or point-like particle systems [20l [29] . In partic- 
ular, we find, at rather large packing fractions, "inverse 
bands" and a possibly direct transition from disorder to 
a "Toner-Tu" [9l4TT[ [25] collectively moving phase. 
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FIG. 1. (color online) (a) Photograph and (b) sketch of 
one polar disk, with particle's polarity n indicated. Typical 
snapshots in the petal-shaped geometry, for the experiment 
(c), and the model (d). Parallel [anti-parallel] particles are 
shown in red (+1) [blue (-1)]. For details refer to [49] . 



The polar disks (Fig. [T^,b) are vibrated between two 
plates. Rather than modeling their full three-dimensional 
dynamics, we describe their effective two-dimensional 
motion. Dictated by the experimental system the main 
new features of the model are: (i) the dynamics of the 
particle's intrinsic polarity with respect to their veloc- 
ity is explicitly described, and (ii) no explicit alignment 
rules are employed, but collisions are modeled explicitly. 
Building on experimental observations, notably that sin- 
gle particles move backward for significant time periods 
with their velocity essentially antiparallel to their direc- 
tor, we were led to the following model: Particle i is 
subject to a noisy acceleration along its polarity axis n* 
(with anisotropic, intrinsic, "active" noise, respecting the 
particle's polar symmetry), balanced by an effective lin- 
ear friction term along its velocity v* = ^r*, with r* de- 
noting the particle's coordinates. Particles i and j with 
I J.* _ J.J I < d^ where d is the particle diameter, interact by 
means of a pairwise, inelastic, repulsive interaction force 
F*-^, yielding the equations: 



d 
dt 



v^ = [/i + 7?||]n^ + T^^ni - /3v^ + ^ F^- 



(1) 



where /i and /3 are constants giving rise to a station- 
ary speed V = /i//3, n^ is a unit vector perpendicular to 
J^\ ^11, ± represent Gaussian distributed white noises with 
zero mean, i.e. (^||,_l(^)^||,±(^0) — 2^||,±^(^ — t'), where 
D\\^j_ denotes the corresponding diffusion constant. The 
interaction force F*-' is given by the established spring 
dash-pot model [50[ [51], which, for hard particles, de- 
pends only a single parameter, the restitution coefficient 

em- 

Eq. (fTl) must be complemented by one governing the 
polarity of particles, which was observed to remain anti- 
aligned to the velocity during episodes of backward mo- 
tion. In other words, when a* = Z(v*,n*), the angle 
between velocity and polarity, is acute, frictional interac- 
tions with the vibrating plate are assumed to rotate n* 
towards v*, while for |Qf*| > 7r/2, n* rotates towards — v*. 
We thus propose the following equation for the polarity 
angle 0* [with n* = (cos0% sin^*)]: 



dt' 



= C sina* sign(cosa*). 



(2) 



where ( characterizes the strength of the coupling be- 
tween polarity and velocity. This parameter is expected 
to be rather small given the observed persistence of n 
even when v changes sign abruptly. 

To make contact with the experimental results, we 
rescale time t -^ t/r^^ with tq the inverse of the vibra- 
tion frequency / = 115 Hz [2 [2^. Length is measured 
in particle diameters (i: x ^ x/(i. Our model possesses 
six parameters: /i, /3, (^ Z^y, D±^ and e. At fixed ex- 
perimental vibration amplitude F, one parameter can be 
eliminated by matching the typical experimental speed 
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FIG. 2. (color online) (a) PDF of ^;|| and (b) PDF of the 
angle a = Z(n(t), Ar(t + r)) [lin-log] for selected values of 
the time increment r. Experimental data are indicated with 
symbols; model data are illustrated with lines. 



with the model's velocity v = /i//3. In the following, we 
use the experimental data gathered at the vibration am- 
plitude F = 2.7, where the most ordered regimes have 
been observed, and for which v = 0.025 jll|2]. 

We first analyze the single-particle dynamics in or- 
der to test the overall quality of the model and to es- 
timate the remaining four parameters (i.e. /3, C, ^||, D±] 
the restitution coefficient e only affects particle interac- 
tions). To find the best-matching set of parameters, we 
consider the following two quantities: the angular dif- 
fusion constant D^j) and the ratio of the displacement 
fluctuations parallel and perpendicular to the polarity 
(definitions see [53 ). Scanning the four dimensional pa- 
rameter space, we select a best-matching parameter set 
for which both quantities agree with the experimental 
value within an accuracy of ±30%. This is approxi- 
mately equal to the imprecision arising due to differ- 
ent preparations of the experimental setup (see [49] for 
more information). In spite of this modest accuracy, the 
model captures quantitatively the observed experimen- 
tal particle dynamics: We compare the distributions of 
the parallel displacements normalized by r, denoted as 
v\\{r) = Ary/r (Ary is defined in [53 ), and of the angle 
a{r) = Z(n(t), r{t + r) — r(t)) to those recorded experi- 
mentally. We find a very good agreement for all values of 
r considered (Fig. [2|. Note that, as expected, the parti- 
cles exhibit backward motion for significant time periods 
(tails in the negative sector in Fig. |2k, and peaks at ±7r 
in Fig.[2]o). 

We now turn to binary collisions, for which the resti- 
tution coefficient e must be chosen. The following re- 
sults are presented for e = 0.4, but we observed that 
changing e in the range ±30% does not influence colli- 
sion properties significantly [49]. Experiments have re- 
vealed that one "encounter" typically involves many suc- 
cessive collisions, where the particles bounce back with- 
out turning their polarity much, so that they quickly 
collide again. These encounters last for a finite time 
and take place over some finite spatial extension. It 
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FIG. 3. (color online) Scatter graph 0[n — Oout for the experi- 
ment (a), and our model (b). Values of the impact parameter 
b are indicated by the color bar. PDF of the duration TcoI (c) 
[lin-log], and the extension icoi of a collision (d) [log- log]. 



was found experimentally that they are well delimited 
using the following criterion: an encounter starts when 
two particles get closer than some threshold collision dis- 
tance, i.e. |r* — r-^l < dc = 1.7, and their polarities point 
"inwards", i.e. |(r^ + n*)-(r^'+n^')| < |r^-r^'| [2^. An en- 
counter ends either when particles are separated by more 
than dc, or their polarities point "outwards". We have 
used the same criterion for our model. Fig. [3] depicts the 
results of a scattering study for the experimental setup 
and the model. Thousands of binary encounters (here- 
after called collisions for simplicity) were recorded, and 
the outgoing relative angle ^out of the two particles plot- 
ted against their incoming relative angle ^in, the impact 
parameter b G [0, 1] [5Ti is shown as color code (Fig. 
[3^,b). The model data shows a striking agreement with 
the results measured in the experiments: most collisions 
actually leave the polarities unchanged (^out — — ^in), 
and a minority of them align the particles almost per- 
fectly (^out — 0)- We estimated the fraction of polar 
aligned events [54J, finding 0.14 for the model and 0.18 
for the experiment. The model also matches the distri- 
bution of head-on (6^0) and glancing (6 ~ 1) collision 
events. We further determined the PDF of the dura- 
tion of collisions TcoI as well as that of their spatial ex- 
tension -^cob given by the center of mass displacement. 
The model reproduces the observed exponential distri- 
bution of Tcoi quantitatively, while it fails to reproduce 
the roughly algebraic decay of ico\ (but nevertheless gives 
a correct mean extension) [55 . 

We performed simulations using the same fiower- 
shaped geometry (Fig. IT]), and number of particles {N = 
890) as in the experiment [ll[2]. For the parameter values 
matching the single particle dynamics and binary colli- 



sions (for vibration amplitude F = 2.7), we observe, as 
in the experiments, fairly large, polar aligned, moving 
clusters (Fig. [l]3,d, for videos refer to [49]). However, the 
order parameter ip{t) = j^\ X^i^ROi ^N? with M{t) de- 
noting the number of particles currently located within 
the central "region of interest" (ROI) of radius 10, is 
typically smaller than in the experiment (Fig. Ilk). The 
effective packing fraction observed in the ROI is found 
to be very close to that of the experiment ((/) c^ 0.39, 
whereas the nominal packing fraction is 0.47), indicating 
that particles accumulate at the boundary in the model 
as well. Running the model at = 0.39 in a box of 
approximately the same size but with periodic boundary 
conditions — a privilege of the in silico approach — yields 
only a marginally larger average polarization (Fig. pk): 
a frustration- free geometry is unable to restore enough 
order. 

We also ran the model in square periodic domains of 
linear size L at the nominal packing fraction (/) = 0.47, 
and then found order being slightly stronger than in 
the experiment (Fig. [4k). Nevertheless, increasing sys- 
tem size L, we observe that the overall order parameter 
(^)^ = (:^l X^i^i ^^\)t decreases first rather slowly, then 
faster (Fig. [IJd, inset). Thus, no true long-range order 
is present at the exact conditions probed experimentally. 
In fact, the correlation length can be estimated by the 
kink in the average polarization as a function of system 
size (inset of Fig.[4]3), leading to a value of approximately 
100, which is larger than the actual experimental system 
size, confirming that order was spanning the whole ex- 
perimental system. 

Next we use a further privilege of in silico investiga- 
tions — the freedom to change parameter values — and 
show that asymptotically ordered regimes would prob- 
ably be observed in slightly different experimental con- 
ditions. Experimentally, the vibration amplitude F was 
used as control parameter for the onset of collective mo- 
tion. Decreasing F to around 2.7 in the experiments, 
order was observed to increase from near-zero to about 
('^)t = 0.5. Unfortunately, due to static friction, the par- 
ticles stopped moving for F values below 2.7. To mimic 
different F- values in the model we multiply both diffusion 
constants D\\ and D± by a coefficient 7^, with 7 G [0, 2], 
so that 7 = 1 corresponds to the experiment at F = 2.7. 
Varying 7, we find the transition to collective motion to 
be close to 7 = 1 (Fig. [4J3). The transition point is ob- 
served to move slightly to the left as the system size is 
increased. This confirms that vibrated polar disks, in the 
experimental conditions, are asymptotically disordered, 
but signals that asymptotically ordered regimes do ex- 
ist nearby, constituting the first report of long-range ori- 
entational order in colliding hard disks without explicit 
alignment. 

Finally, we have performed a systematic exploration of 
the model varying 7 and the packing fraction (j) in square 
domains of linear size L = 200 with periodic boundary 



conditions (Fig. ^). For ^ < 0.6, varying 7, we observe 
the usual phenomenology of models with (effective) po- 
lar alignment like the Vicsek model ffl [T^ [TSl [2Ql [^: 
immediately below the transition, the particles spon- 
taneously segregate in high-density high-order "bands" 
traveling in a low-density disordered sea (Fig. Hrl). Fur- 
ther away from the transition, these nonlinear structures 
disappear, leaving a statistically-homogeneous Toner- Tu 
phase with its characteristic giant number fluctuations 
and long-range correlations [9 HTT| |25 ] . However, we de- 
tected, for large enough packing fractions, narrow disor- 
dered channels (see Fig. [4]i, (4)) for small noise values 
(green circles in Fig. pfc). These "inverse bands", not 
found in dilute or point-like particle models, seem to co- 
exist with the Toner- Tu phase. We believe that the in- 
creased frequency of collisions at large packing fractions 
trigger the emergence of these inhomogeneous structures. 

Interestingly, for (j) > 0.6 we could not observe bands 
(Fig.jIJ^). This suggests a possible direct transition from 
the disordered to the Toner- Tu phase. At this stage, 
however we cannot conclude, due to numerical limita- 
tions, whether this feature remains in the limit of large 
system sizes and asymptotically large times: the width 
of the bands increases with increasing ^ (cf. Fig. [4]i) so 
that their disappearance might just be a finite-size effect. 
However, the longitudinal density profile around (j) ~ 0.6 
turns out to be rather flat, with an overall rather low 
order (as low as {ilj)t ^ 0.2 for (j) = 0.6 and 7 = 1.4). 
They may thus be of different nature from the Vicsek- 
like, sharp, well-ordered bands found at low 0, and could 
cease to exist asymptotically at a packing fraction below 
the rise of jamming and crystallization effects. 

To summarize, we have built a simple yet quantita- 
tively faithful model for the dynamics of the vibrated 
polar disks studied in [U |2] . This model constitutes one 
of the first in which the dynamics of the particle's intrin- 
sic polarity with respect to their velocity is taken into 
account [56, 57 . An adequate description of the granular 
system of vibrated discs requires accounting for the po- 
larity as a slow variable compared to the velocity, which 
can change fast due collisions with the plate or neighbor- 
ing particles. Our in silico study has shown that in the 
original experiments the most ordered state reached was 
in fact in the region of the transition to collective motion, 
slightly on the disordered side. However, asymptotically- 
ordered regimes do exist nearby. The new features of the 
phase diagram, i.e. the emergence of "inverse bands" in 
the low noise regimes of sufficiently dense systems and 
the possibility of a direct transition from disorder to a 
collectively- moving Toner- Tu-like phase, deserve further 
investigations. In particular, this last point, if confirmed 
in the future, might reopen the debate about the possi- 
bility of a continuous transition to collective motion since 
the structures "responsible" for its discontinuous charac- 
ter — the bands — would then not exist. 
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FIG. 4. (color online) (a) PDF of the average polarization 
ip, evaluated within the ROI, for the experimental system, the 
model in the petal-shaped geometry and in periodic bound- 
aries using two values of packing fractions: = {0.39,0.47}. 
(b) Average polarization {^)t as a function of the noise frac- 
tion 7^ = i^||/i^|j^^^"^ = D^/D^^^"^, shown for three boundary 
sizes L e {50,100,200} and = 0.47. Inset: (^)t [log-log] 
for 7=1 and = 0.47 as function of system size L. (c) 
Sketch of packing fraction(0)-noise(7) phase diagram: States 
with {^)t < 0.5 are indicated by ■, polar homogenous states 
with (^)t > 0.5 by A, and states exhibiting heterogenous pat- 
terns transversal to the average moving direction ("bands") 
are depicted by •. (d) Representative snapshots for selected 
0-7" values indicated by numbers in (c). 
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^J — \ri-r3\ ' 

and parameters A and k are related to the (normal) 
restitution coefficient via e = exp (— Arhc) with The = 

^0 = ^{[Ht + ^) - Ht)?)t (with (...)t denoting a 
time-average), and the ratio of the displacement fluctu- 
ations parallel and perpendicular to the polarity, H := 

^^friyl , with Aril x(r,i) = E^'^lV^/lAmAt + 

At)-r'(mAt)]-n||^^(mAt), where At = 5.75ro, equal to 
the experimental sampling time determined by the cam- 
era frame rate. 

To avoid counting events whose outcomes are possibly 
influenced by the limited camera frame rate, we neglected 
all events with Oin < 7r/6. Moreover, polar events are 
defined by: ^out > -^in/2. 

To what degree this is an actual discrepancy between 
model and experiment remains to be clarified. In fact, 
the very existence of an algebraic decay for the experi- 
mental data can be questioned due to the small number 

of collisions with large extensions. 
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